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Abstract 


Skidding is one of the most dangerous instabilities which a vehicle might encounter 
while braking. This phenomenon becomes all the more critical in adverse conditions 
like flooded roads, roads covered with sheet ice etc. The importance of anti-skid land- 
ing gear is also highlighted by the stringent operational requirements of the defence 
forces. 

In the present work a simple anti-skid system has been presented. The equations 
are for a general case of the vehicle problem but the numerical solution has been 
carried out for a single wheel hcavc-degrecs-of-freedom vehicle model on a linear 
spring and linear damper shock absorber. The dynamics of the vehicle during anti- 
skid braking have been presented. 

The anti-skid system has been implemented by taking into consideration the 
stochastic properties of the road. The braking force is predicted one step ahead 
in time with the help of an auto-regressive or an all poles filter. The road is mod- 
elled as the combination of a deterministic variable mean process and a random zero 
mean process. The road is modelled in the space domain, and transformed into the 
time-domain with the help of an integral space-time equation. The road is simulated 
in terms of the sum of cosine functions with random frequencies and random phase 
angles. 
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Chapter 1 


Introduction and Literature Survey 


The present work deals with the optimal braking of aircraft landing gears during 
the landing run. An attempt has been made to achieve some amount of anti-skid 
behaviour in the braking system by using some basic tenets of control theory and 
theory of stochastic processes. The salient aspects of this problem, review of the 
available literature, and the approach adopted in the present study, follow in the 
subsequent sections of this chapter. 

1.1 Introduction 

In the modern world, the transportation sector, and therefore, vehicles play a major 
role in the day-to-day life. The two principal requirements from a transportation 
system, other than its obvious requirement to carry goods or passengers, are the 
abili ty to accelerate/decelcrate at will, and the ability to change its direction of motion 
at will. The acceleration is provided by complex propulsion mechanisms, where as 
the deceleration is carried out by relatively simplistic mechanisms. 

The research on driving mechanisms with its obvious economic implications has 
been extensive, however, that on the braking mechanisms has languished to some 



extent. Though this is true to a certain extent for land vehicles, for aircrafts, the 
braking mechanisms have indirect economic implications due to their ability to dictate 
operational requirements pertaining to landing. Stark examples of this are the present 
day “jumbo jets ”, which require large and well prepared runway surfaces and are 
therefore restricted to operating between major airports. It can be easily seen that 
the present trend of growth in air traffic will lead, in the future, to requirements of 
larger aircrafts landing on smaller landing strips. 

In the case of the military, the requirements from a landing system are much more 
severe and diverse. The military aircrafts often have to takeoff and land on badly 
prepared small runways. This necessitates extremely sturdy and specialised landing 
gears and brakes. On aircraft carriers the landing deck is often slippery due to the 
shipping of green seas in rough weather, and the runway length quite limited. This 
necessitates specialised mechanisms to effectively halt the aircraft. 

For land vehicles, the design of a braking system is often dictated by stringent 
safety requirements which require a minimum braking distance and time. However as 
can be seen in §1.2 the braking distance can not be reduced arbitrarily, as doing so 
would have implications on the general dynamics of the vehicle and it’s safety. The 
main danger, specially on bad roads like those covered with sheet-ice, snow, water, 
etc., is skidding. 

The above scenario points to the needs of an u intelligent ” braking gear capable 
of reducing braking distances drastically on any kind of surface. 


1.2 Dynamics of Braking 

1.2.1 The Limitation on Maximum Braking Force 

In a wheeled vehicle the process of braking involves a change in the dynamics of the 
moving body which if not checked could proceed into unstable regions. 
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Any rolling body like a wheel has inherent stability due to the stabilising gyro- 
scopic forces set up due to its motion. 1 hese forces tend to oppose the effects of any 
destabilising moments that may exist due to unbalanced or misaligned forces on the 
body. 

The rolling motion of the free rolling wheel is guaranteed by the frictional force 
between the wheel and the track. Thus if a braking force greater than the maximum 
frictional force sustainable by the road is applied on the wheel, it stops rolling and 
starts skidding on the surface. In the case of total skid the stabilising gyroscopic 
forces are lost and the destabilising forces come in full play. This leads to a loss of 
control on the vehicle and may eventually lead to an accident. Further, any corrective 
measures other than the release of brakes may be to no avail, as they may not be 
sufficient in magnitude to overcome the destabilising moments on the vehicle. 

1.2.2 Variation of Maximum Friction Force with Time 

As the vehicle moves over the track the effect, of its weight passed on to the track, 
filtered through the landing gear, varies due to the dynamics of the vehicle. The vehi- 
cle oscillations induced by the random unevenness depend on the vehicle suspension 
characteristics, its forward velocity and the nature of the track profile. The temporal 
variation in the filtering of the weight through the landing gear naturally results in 
variation of the reaction force by the track on the wheel, and a corresponding vari- 
ation in the maximum sustainable friction force. Both these forces are random in 
nature due to the randomness of the track profile. 

1.2.3 Variation of Frictional Coefficient 

The frictional coefficient for a rolling tyre and the road is not uniform but varies as 
in Figure 1.1 [1], with the skid. It has a maxima somewhere between 0.05 and 0.3s 


v 
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Figure 1.1: Variation of friction coefficient /x with tyre skid 


where slip s is defined as [1] : 


s = 1 


where u; is the angular velocity of the wheel and u a the free roll angular velocity 
at that speed. This implies that in order to have better braking characteristics, the 
vehicle has to continuously vary either the braking force or the vehicle parameters 
so as to change the filtering characteristics for the better, and maintain itself in the 
range of optimal skid. 

\ 

1.3 The Control Problem 


The generalised vehicle control problem may be summed up as the modelling of the 
vehicle, its input disturbances or forces (track parameters, the drag and driving forces) 
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and its response to these forces. Control on the response is effected by varying either 
the model parameters or the force parameters, whichever may be controllable. It 
is also possible to optimise for the model parameters for some suitable optimality 
criterion. The optimisation (say of the mean square response) can be carried out for 
a statistical description (spatial or temporal) of the disturbance parameters. This 
will give an optimised response for a probable range of disturbances the vehicle might 
encounter in reality. Such optimisations lead to better and more economic designs of 
the vehicles. 


1.4 Literature Survey 

As described in §1.3 the control problem of a vehicle deals with 3 major modelling 
tasks namely 

• modelling of the input process, 

• modelling of the system, and 

• the control mechanism 

Extensive work has been done on each of the above, and considerable literature is 
therefore available, some of which is detailed below. However literature on anti-skid 
braking was not discovered. Two such systems are commercially available, developed 
by Messrs. Goodyear and Dunlop. This problem has also been studied by various 
defence establishments. The lack of literature is perhaps due to the strategic nature 
of this work. 

1.4.1 Input Models 

The inputs in a vehicle problem are chiefly the track undulations, therefore the mod- 
elling of the input is modelling of the track. The track is usually modelled as a 
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discrete or continuous random process. Initially the track profile was taken directly 
from observations on roads, runways or rail tracks. Tung et al. [32] have measured un- 


dulations of runway profiles for their input models. The track can also be completely 
described by its statistical properties like the auto- correlation and power spectral 
density (PSD) functions to the second degree. Lindgren [20] cites that the PSD can 
be given as 



Sr<Z et 


0 


(/ ^’min ^ W i ^rnax 

otherwise 


( 1 . 2 ) 


where u) is the temporal frequency, w rr j = 2n, u and u max are the cut-off fre- 
quencies. He has modelled the track as an auto-regressive process of order 3 given 


as 


q(s + h ) 

sf(“0 


aiq(s) + a 2 q{s - 2 h) + a 3 q(s - 2 h) + e(s) 

( ha l\ 1 


2ir 


|1 — ci\e^^ h — a 2 e“ 2jfu;/l 


a 3 e 


—3jujh [2 


(1.3) 

(1.4) 


Yadav and Nigam [35] have modelled the track profile as a homogeneous process with 
zero mean and the single sided spatial PSD function 

U2+1/2 */ D < Do 

0 otherwise 


<*>(0) 


(1.5) 


cf}{w) 


where Do is the cut-off frequency and v is the correlation constant. The track was 
then generated as the output of a cosine filter to white noise. ElMadany [6] has used 
a single sided temporal PSD given as 

' {—■) •/" < «• (i 6) 

0 otherwise 

where v is the vehicle forward speed, a 2 is the variance and or a coefficient depending 
on the road irregularities. The excitation model was then obtained by a shaping filter 
to white noise 

A(<) = -ovX(t) + £(i) • (1.7) 
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Narayanan and Raju [21] have described their track by the auto-correlation func- 
tion 

R(0 = <7 2 e~™ (i. 8 ) 

where — $2 — S\ : is the spatial lag with a > 0. The effect of the rolling contact of 
the tyre has been considered as an additional filter in cascade [15]. 

h"(s) + (o -f jS)h'(s) + a/3h(s) = kiv(s) (1-9) 

expressed in the time domain this is, 

h\ = -{a + 0)s{t)h\ - a/3s{t)hi + -i(<)w{a(<)} (1.10) 

a 

where o is the cut-off wave number of the road profile spectrum, (3 is the cut-off 
wave number of the rolling contact filter, k = i3cr^2a is a factor to ensure the 
desired variance, and cu{s(t)} is a white noise process transformed in time domain. 
Freymann [10] in his experimental investigation has described the runway as a series 
of bumps and has devised a bump factor 

= i/T+1| + 2eIc2 (1.11) 

where E& = e - ^, 

Ca = coscu^a/I - £ 2 , and 
= cooA/V 

where A is the bump spacing, £ the damping ratio, u 0 the natural frequency of the 
landing gear, and V the vehicles forward speed. He has gone on to investigate the 
most critical bump factor. 

1.4.2 Vehicle Models and Solution Schemes 

The increasing demand for safety and ride comfort, especially at high speeds, has led 
to the development of new types of suspension systems for road and rail vehicles. 
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Sharp and Crolla [25] have reviewed different concepts in road vehicle suspension 
system design, and have compared relative merits and demerits of passive, active, 
semi -active, and slow- active suspensions from the view point of ride comfort and 
control of wheel load variations* In recent years a lot of research work has been car- 
ried out in the analysis and design of active suspension systems. Tomizuka [30] has 
studied the preview control of vehicle suspension from the viewpoint of discrete opti- 
mal control. II ac [11, 12] has used stochastic optimal control theory in the suspension 

a 

optimisation of a two -degree-of- -freedom vehicle model and a vehicle with an elastic 
body, traversing a rough road, for better ride comfort and road holding characteris- 
tics. Yoshimura et al. [36] have considered a cascade arrangement of a Kalman filter 
and optimal controller, to formulate and solve the problem of designing an active sus- 
pension to control the vertical vibrations of railway track vehicle. He has performed 
measurements in a noisy environment and estimated unknown state variables from 
this data by using a Kalman filter. Yoshimura and Sugimoto [38] have continued the 
same work for a vehicle travelling on a flexible beam wit h an irregular surface. The 
vehicle is modelled as a four degrees of freedom system. The beam defection itself 
is described by a set of modal functions. The measurement of track inputs is done by 
means of a pre-view arrangement. ElMadany [6] used stochastic optimal control and 
estimation theories to design active suspension system for a cab ride in a tractor- 
semitrailer vehicle. Sharp and llassan (26] have shown that a pneumatic air-pump 
based automobile suspension system has very good ride performance properties over 
a tvide range of workspace road roughness ratios. Lindgren [20] has analysed the 
response of a single- degree of freedom system for the case of vehicle jump or take-off 
while traversing a rough road. lie has presented a switching stochastic differential 
equation for the system. Thompson [29] has optimised the tuning and damping ratio 
for a two-degrees-of-frecdom system, to minimise wheel bounce on a random road. 
Freymann [10] has carried out experimental investigations for aircraft taxi ground 
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loads alleviation of landing gear. In these studies, the vehicle was assumed to be 
traversing the rough road/track with a constant velocity. The stationary response 
of the vehicle was found to be superior in comparison to that of a passive system in 
terms of ride comfort and road-holding. 

Non-stationarity is inherent to most vehicle problems. The vehicle response is 
non-stationary 

• for variable velocity traverse over rough roads with spatial homogeneity, 

• when the vehicle traversing a smooth surface (zero initial conditions) at con- 
stant velocity suddenly encounters a spatially homogeneous rough terrain and 
continues to move — in this case the response becomes stationary after a lapse 
of time, 

• for traverse over a partially noil-homogeneous terrain with constant or variable 
velocity. 

Non-stationary vehicle response due to variable velocity has been treated in a 
time-domain formulation by Virchisand Robson [33] and by Sobczyk and Macvean [28] 
A space-domain formulation admitting a evolutionary spectral form in the wavenum- 
ber has been presented by Yadav and Nigam [35]. In their paper the knowledge of 
the exact nature of the time space relation was a pre-requisite. They transformed 
their temporal equations into the space domain assuming the time space relation to 
conform to a power law. Hammond and Harrison [13, 15] adopted linear state-space 
techniques, yielding a set of matrix differential equations for the determination of the 
non-stationary response covariance. In the works mentioned above, a single-degree- 
of-freedom vehicle model with passive suspension was assumed, and the different 
analyses involved evaluation of multiple integrals in the time-domain, contour in- 
tegration in the wavenumber domain, and numerical integration for the covariance 
response. 


9 



1.5 Present Work 


The present work is an attempt to design an anti-skid algorithm for an optimal 
braking problem. 

In aircraft braking systems during landing there already exist two types of anti- 
skid methods. 

1. Pulsed braking system : 

In this system the braking force is applied in discrete pulses. The wheel is 
braked for a short period of time with no concessions to skidding, and then 
released. During the period of free roll the wheel recovers from any skid it 
might have entered. This process is repeated until the aircraft is brought to a 
halt. This is a rugged design even though very sub-optimal. 

2. Skid monitoring system : 

In this type of landing gear the motion of the aircraft is monitored to detect 
any skidding during braking, with the help of sensors. When a skid is sensed, 
the wheel is released to roll freely and come out of the skid. This process is 
repeated until the aircraft is brought to a halt. 

In both the above methods the braking force is either switched on or off but not 
varied continuously. The present work attempts to optimise the braking by varying 
the braking force with time. 

It is physically impossible for any device to react instantaneously to a stimulus. 
This is mainly due to the time taken : 

• to detect the stimulus and its nature. 

• to calculate the response. 
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• by the device to activate it’s sensors and achieve fully the response, due to the 
inherent physical inertia of the servo-mechanism. 

Thus a continually varying braking force will lag behind the impulse by a small 
time interval. Therefore to achieve optimality it is not sufficient to just find out 
the braking force and apply it, but to do it with a view to the future. That is, 
the brakes have to be applied so that, the events, namely the vehicle encountering 
the impulse and the brakes achieving the calculated optimal braking force, should be 
simultaneous. This requires the prediction of the impulse one step ahead in time. The 
prediction of the impulse can be done by using a pre-view arrangement as suggested 
by Yoshimura [37, 38] or by using statistical and/or probabilistic methods. 

In this study the surface is assumed to be a random process with a given PSD. The 
surface characteristics are predicted one step ahead in time by using the statistical 
properties of random processes, with the help of an auto regressive filter. 

1.5.1 Mathematical Models Used 

This study is computational in nature and deals with mathematical modelling of the 
vehicle and the track surface. 

The Vehicle Model 

The vehicle is modelled as a rigid body, single wheel, heave degrees of freedom model, 
supported on a linear shock absorber consisting of a spring and a damper. The tyre 
is modelled as an equivalent linear spring and damper as in Figure 1.2 

The Surface Model 

The surface is modelled as a zero mean random process superimposed over a variable 
deterministic mean. The PSD of the zero mean random process has been found to 
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Figure 1.2: The vehicle model 


conform to a “flat bell shape” . Various approximations to this shape have been used 
to describe the PSD distribution, depending on the primary variable chosen, namely 
velocity or distance. In this work the PSD is described by equation 2.5 [22] : 


*(«) = 



0 


if jO j < D 0 

otherwise 


Here $ is the PSD, D the spatial frequency, Do the cut-off frequency, a is the standard 
deviation of the sample track from the zero mean, and oj is a correlation constant 


describing the roughness of the track. 

The random process is generated using a cosine filter [27] to a uniform deviate, 
described in §4.1. The uniform deviate is generated using the modulo, and the data 
encryption standard (DES) algorithm [24]. The track heights generated from this, 
as a function of distance, are added to the corresponding deterministic means which 
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were taken as simple algebraic or transcendental functions of distance. 

Prediction Algorithm 

The auto-regressive algorithm was used for predicting the track profile one step ahead 
in time. In this algorithm the predicted value is taken to be a weighted average of 
the time series’ history, as in equation 4.9 [34]. 

< = a kXn-k (1-12) 

*=i 

Where x* is the predicted value at the present instant, a* the weighting factor, and 
x n -k the actual value of the time series at an instant k time steps in the past. Here 
ajt are calculated such that the mean square error in prediction is a minimum. Er- 
ror is defined as the difference between the actual value of the time series, and it’s 
corresponding predicted value at any instant of time. 

The equations of motion are solved in the time-domain using the impulse response 
method. This method calculates the response at the present time due to a unit impulse 
disturbance in the past. This response is convolved with the disturbance history to 
arrive at the present response of the system. 

The programmes are developed in FORTRAN-77 on the IIP— 9000 series 850 com- 
puter systems with a UNIX operating system, at I. I. T. Kanpur. 


i 
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Chapter 2 


Models Chosen 


The vehicle considered for the present work is a light fighter / trainer. Only a single 
wheel landing gear has been considered. The models for the aircraft and the runway 
profile chosen for the present work are detailed in the subsequent sections of this 
chapter. 


2.1 The Vehicle Model 

The vehicle is modelled as a rigid body, single wheel, heave degrees of freedom system, 
suspended on a linear spring and damper shock absorber. The weight of the landing 
gear is considered to be partially sprung and partially unsprung. The sprung mass of 
the landing gear is considered along with the vehicle mass. The wheel and the tyre 
are modelled as an unsprung mass. The damping and spring effects of the tyre are 
modelled as an equivalent linear spring and damper. The model is as presented in 
Figure 1.2. 



2*1.1 Forces on The Vehicle 


The forces acting on the system which are considered in this work, other then the 
track undulations, are 

1. Lift force on the aircraft, 

2. Drag on the aircraft, 

3. Gravitational force, in the form of weights of the sprung and the unsprung 
masses of the aircraft and landing gear, 

4. Braking force being applied on the aircraft, and 

5. Reaction of the ground as visible to the tyre. 

Gyroscopic moments/ forces on the wheel are neglected, as they do not have a direct 
bearing on this problem. Loss of energy in the shock absorber due to heating, sec- 
ondary vibrations such as; vibrations of the shock absorber wall, emission of sound 
etc., are assumed not to effect the damping and spring behaviour of the shock ab- 
sorber, and are neglected. Loss of energy due to bearing friction and other such 
secondary phenomenon is also neglected. The system is assumed to be isolated and 
adiabatic. ' 


2.2 The Input Process 

The input disturbance to the vehicle model is the track unevenness, which gets fil- 
tered successively through the tyre and the shock absorber. The track unevenness is 
modelled as a random process in the space domain, as described below. 

The track is assumed to consist of two distinct processes namely, 

• a deterministic variable mean process, and 
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• a random zero mean process described by it’s statistical properties. 

The track height is the algebraic sum of the two sub- processes. 

2.2.1 Tlie Deterministic Variable Mean Process 

This process can be described by infinitely many combinations, over distance, of 
simple algebraic or transcendental functions, to satisfy some physical properties of a 
given runway like 

• differential settlement of concrete slabs, 

• bomb craters covered with repair mats, 

• inclined runways, 

• parabolic or curved runway surfaces, 

• unevenness in landing decks of aircraft carriers due to buckling, etc. 

In this study three cases of the runway profile have been considered namely 

• a stepped runway, namely a crater covered with a repair mat, 

• an inclined runway, and 

• a sinusoidal runway. 

2.2.2 The Random Zero Mean Process 

This process can be described either by it’s auto-correlation function or the Power 
Spectral Density [PSD] function. 
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In the literature surveyed most of the work has been done using the following 
approximation for the auto-correlation function for the road or it’s corresponding 
PSD function [22]. 


R{t) = <rV°M (2.1) 

m = ~~~T7 — 2 ( 2 - 2 ) 

7T W l -fa 2 

However as can be easily seen, the auto-correlation function is not differentiable at 
t = 0, therefore X(t), the corresponding random process, is also not differentiable. 
The following alternative approximation function for the auto-correlation, however is 
differentiable for all r [22]. 


2 ^ — OfT 2 


R(r) = <r 2 e 

The corresponding PSD function for the auto-correlation equation is 

2 


<!>{“) 


c r 


2 v / 7ra 


e 


(2.3) 


(2.4) 


Narayanan and Raju [21] pointed out that the road may be treated as a homoge- 
neous random process only in the space domain. Therefore in this study the runway 
has been generated in the space domain and then converted to the time domain via 
a space-time relation. For generating the random process therefore the spatial PSD 
and the spatial frequency have been used. Also since a discrete process is generated, 
the sampling theorem dictates a cut-off frequency, within which region only the PSD 
need be integrated. Thus equation 2.4 becomes 


*(«) 


T^ C 1 


(2.5) 


if |S![ < !!„ 

0 otherwise 

where Cl 0 is thd cut-off frequency given by [24] 

n » = 4 

where A is the sampling interval. In this study however, D 0 was chosen so that $(fl) 
was sufficiently small. A was calculated for this value of 0 0 - 


( 2 . 6 ) 
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Chapter 3 


Equations of Motion 


The equations of motion for the model described in the previous chapter, forces acting 
on it, and the physical restrictions imposed by the shock absorber and the tyre, are 
detailed in the subsequent sections of this chapter. A detailed description of the 
changes in the model due to the tire/shock absorber locking during their operation 
is also presented. 


3.1 The Generalised Equation of Motion 

The dynamics of a rigid body can be described adequately by the following system 
of ordinary differential equations. 

MZ + CZ + KZ = F (3.1) 

Where M is the inertia matrix, C is the damping matrix, and K is the stiffness 
matrix. F is the vector of external forcing functions, and Z the state vector or the 
displacement vector. 



3.2 Solution for The Equation of Motion 


The solution of equat ion 3. 1 is got by using t he impulse response method [5, 19, 22, 31]. 


Defining a vector Y = 


Z 

z 


equation 3.1 can be rewritten as 


Y + AY = P 


(3.2) 


Where 


A = 


M _I C M -1 K 

-I 0 


L 0 J 

Where I is the identity matrix, and 0 is the zero matrix/ vector. 

Equation 3.2 is in the standard form of the Bernaulli equation [17]. Thus it’s 
solution can be written as [17] : 


Yef Adl — J Pef Adt dr -t- C, 


(3.3) 


Here Y, A, and P can be time dependent and Ci is a constant depending on the 
initial conditions. However for a linear system A is a constant with time. Thus for a 
linear system, within the integral limits of t 0 and t, equation 3.3 becomes : 


[ A dt [‘ , , f Adi j , _ 

Ye J, « = / P(r)e J 'o dr + Ci 
Jt 0 


y e A-(i-i 0 ) = J* p( r ) e A -( T -‘°) dr + Ci 

Y = [ l P(r)e A (r - f) dr + Ci- e " A (< - to) 

Jto 


At t = io equation 3.4 gives 



(3.4) 


(3.5) 
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Substituting equation 3.5 in 3.4 and defining the impulse response function 

h (f) = e~ A< ' (3.6) 

We can write 

Y(i) = Y(/ 0 )h(Z - t Q ) + f h (t - r)P(r) dr (3.7) 

J to 

The evaluation of the exponential in equation 3.6 can be done as follows : A can be 
expressed as 

A = UEtr 1 (3.8) 

where U is the modal matrix of A, and E is a diagonal matrix whose diagonal elements 
correspond to the eigen values of A. 

It is easily seen from 3.8 that 

A 2 = UEU -1 UEU _1 
= UE 2 U _1 

A 3 = UE 2 U“ 1 UEU" 1 
= UE 3 U _1 
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From 3.10 and 3.6 we get 


h(<) = UVU - 1 


Where U is the modal matrix of A, and 


V = 



V i = j 

V i ± j 


(3.11) 


where aj are the eigen values of A. 


3.3 Present Model 


For the present model the matrices M, C, K, F, and Z are as follows 


M = 


m i 0 

0 77 ? 2 


c = 


Cl -Cl 
-C, Ci + C 2 


i 


K = 


K x -K x 
-K x I<i+I<2 


F = 


m 

C 2 h + I< 2 h 



The sign conventions are as shown in Figure 1.2. Z\, Z 2 , h, and t are assumed to 
be zero at touch down. 
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3.4 Forces on The Model 

As enumerated in §2.1.1 the external forces acting on the aircraft are Lift, Drag, 
Gravity, Reaction, and Braking. The action of gravity is incorporated in the general 
equation via the mass of the vehicle. The expressions for the other forces for the free 
rolling case are as below. 


C L = 


(3.13) 


3.4.1 Lift Force on The Aircraft 

The Lift on the aircraft can be represented by the following cardinal form [3] : 

L = l -pV 2 SC L (3.12) 

where p is the density of air, V is the speed of the aircraft, S is the lifting surface 
area or the area of the wing, and Cl is the co-efficient of lift given by 

Cl 0 + &Cl a V a < Q stall 
0 V a > Qstall 

where a is the angle of attack, Ostiii is the angle of stall, Cl 0 is the lift coefficient at 
zero angle of attack, and Ci a is the lift curve slope. For most aircraft wings the lift 
varies linearly with the angle of attack, and so the lift curve slope is a constant for 
the operation range. The angle of attack varies as the aircraft taxiis due to the heave 
of the wings. The effective angle of attack is got as 

-i (Zi\ 

Q c/J — °o + tan I — I 

where c*o is the pre-set angle of attack of the wing. 


(3.14) 


3.4.2 Drag Force on The Aircraft 

The Drag on the aircraft is also expressed in a cardinal form similar to 3.12 as follows 
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where Cd is the drag coefiicient, which is usually described by the “drag polar” below 

Cd = Cd 0 + kCl (3.16) 

The total drag is got by summing up the drags of individual elements like the wing, 
fuselage, tail plane, etc. Lach of these individual drag forces can be represented in 
the form of equation 3.15. 

3.4.3 Reaction of The Ground 

The weight of the aircraft is filtered through the tyre, and the landing gear. The 
weight as visible to the track is 

R = C 2 (h - Z 2 ) + K 2 (h - Z 2 ) (3.17) 

This weight manifests itself in the form of the reaction of the ground on the vehicle. 


3.4.4 Braking Force Applied on The Aircraft 

The braking force on the aircraft is limited by the dynamics of the vehicle as shown 
in §1.2.1. The maximum applicable braking force can be given as 

B = // max 7? (3.18) 


where fi is given as a function of skid s, and has a maxima between 0.05 and 0.3 as 
shown in Figure 1.1 [lj. Skid s is defined as in equation 1.1 


1 


u; 

k-’a 


Where u> is the angular velocity of the wheel and 

V V 


U>a = 


rolling radius ~ — - (3.19) 

The friction also depends on the vehicle speed, tyre inflation pressure, tread depth 
and ground surface conditions [1]. 



3.5 Domain of The Equations of Motion 

The generalised equations of motion, and the corresponding equation for the reaction 
force are not valid for the following cases : 

• tyre bottomed, that is — Tj < Z 2 — h < T e , 

• shock absorber locked, that is —S\ < Z\ — Z 2 < S e , 

• the aircraft ballooning or jumping/bumping, that is R < 0. 

3.5.1 Equations for Tyre Bottoming 

Tyre bottoming occurs when the rim starts touching the ground. Even though it 
might be allowable for an instant, on larger periods this may result in the tyre burst- 
ing. Tyre extending is normally accompanied with the ballooning of the aircraft. 
Pure extension of the tyre without the aircraft is a hypothetical borderline case. 

In the case of tyre locking the tyre looses it’s filtering properties and the track 
disturbances are passed directly to the shock absorber. The system now behaves like 
a one-degree-of-freedom -system riding on a mass which moves along the ground. 
The equation becomes : 

mi Z i + Cx(Z x - Zi) + I<i{Z] - Z 2 ) = L - m ig 

or 

rriiZi + C\Z\ + K\Zi — L — rn\g + Ky Z 2 -f C\Z 2 (3.20) 

here due to the lack of tyre filtering Z 2 — h + 6 , Z 2 = A, and Z 2 — h. 6 is the 
displacement of the tyre at which it is locked. 

The reaction on the ground is now given by the following equation : 

R = m 2 (Z 2 + g) + A'l (Z 2 — Zi) + C\{Z 2 — Z\) + Flocking (3.21) 

where Faking is the force required to lock the tyre. 
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• 5.2 Equations for Shock Absorber Locking 

he locking of the shock absorber is not as critical as the locking of the tyre. In 
iis case the filtering characteristics of the shock absorber fail and the aircraft starts 
ding on the tyre, so that the tyre accelerations and displacements get transmitted to 
, directly. This results in very high instantaneous accelerations of the aircraft, which 
ray be uncomfortable to the passengers. If the aircraft does not recover quickly, the 
ituation may lead to the locking of the tyre signifying, a total collapse of the landing 
ear. 

The equations of motion for this case are : 

m 2 Z 2 + C 2 (Z 2 — A) + K 2 (Z 2 — h ) = rriiZi + (mi — m 2 )g + Fi oc ki n3 ~ L 


>r 

m 2 Z 2 -f- C 2 Z 2 + I< 2 Z 2 = mi'Z\ + (7721 — m 2 )g + C 2 h + K 2 h + Fi oc ki n g — L (3.22) 

lere due to the failure of the shock absorber Z\ = Z 2 + 5, Z\ — Z 2 , and Z x = Z 2 . 6 
s the displacement of the shock absorber at which it is is locked. 

The equation for the reaction force does not change as it depends only on the tyre 
parameters. 1 

3.5.3 Equation for Aircraft Ballooning 

The ballooning of the aircraft takes place when the runway is very rough and dips 
sharply after a peak or in the early stages when the aircraft is recovering from the 
impact of touch-down. The centre of gravity of the aircraft — landing-gear system 
follows the trajectory of a simple projectile, however the wheel and the aircraft con- 

I 

tinue to oscillate with respect to one other. The governing equation is that of the 
two-degree-of-freedom system with the external forcing function zero. The reaction 
on the ground in this case is obviously zero. 
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3.6 Space— Time Relationship 


In this study the track has been generated in the space domain, where as the brak- 
ing and the dynamics are being studied in the time domain. This necessitates the 
determination of the space-time relationship of the system. The relationship is not 
easily definable mathematically as the braking force is continually varying with time 
in a random fashion. However it can be stated as an integral equation which can be 
integrated numerically in real-time. The relationship can be stated as follows : 


s(t) = 


D + B 


mi + m 2 


r*- D + B 


' t 0 mi + m 2 

rt 


dt + s(to) 


m = - / 

J k 

s(t) — — f s(t) dt + s(t 0 ) 

J to 


where D is the drag force and B is the brake force. 


(3.23) 

(3.24) 

(3.25) 


3.7 Skid— Time Relationship 


The skid is defined in equation 1.1, as a function of the wheel’s angular velocity. 
The angular velocity in turn depends on the variation of the braking torque, and the 
moment generated due to the shift in the centre of pressure of the radial load. The 
angular acceleration of tfie wheel is thus given as : 






(3.26) 


where D is the free outer diameter of the tyre, 8 is the tyre deflection, X c is the 
forward shift in centre of pressure, and Tj is the braking torque being applied. 

The angular velocity can thus be got by integrating equation 3.26 in the time 
domain. 

w — f dj^t') dt + oj(to) (3.27) 
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Chapter 4 


The Track Profile Generation and 

Prediction 


In this chapter the method used to generate the track profile, described in chapter 2, 
is detailed. The prediction routine for forecasting the track height, one time step in 
the future, is also discussed in the subsequent sections of this chapter. 

4.1 Generation of The Track Profile 

The track profile is described in §2.2. The deterministic variable mean process is easily 
generated as it is. fully described by it's equation, whereas the zero mean random 
process has to be generated using it’s statistical properties. The random process 
remains homogeneous in the space-domain and becomes noil-homogeneous in the 
time-domain. Thus it’s simulation in the space-domain is simpler. 

The random process is simulated in terms of the sum of cosine functions with 
random frequencies and random phase angles [27] as below 



where 


(4.2) 


fOO I J 

a ~ / 4>o(u>) dui 

is the standard deviation of the process f 0 (t), u k (k = 1,2 are independent 
random variables identically distributed with the density function g{u>) = g{u>k) ob- 
tained by normalising ^o(tu), 


(4-3) 

cr* 

and tpk are independent random variables identically distributed with the uniform 
density between 0 and 2 tt. tck and ipi(k,l = 1,2, • • • ,-/V) are independent. 

The random process t hat is generated can be easily seen to be of zero mean and 
the required PSD [27]. The ensemble average £'[/(<)] of f(t ) is 


mo] 


a 


($)’ r2* 

/ / cos(cet + il’)g(bj) dxp du> = 0 (4-4) 

2 7T J-oo Jo 


The auto correlation function R(t) of f(t) is 


R{r) -■ E[f(t -f r ■)/(/)] 

‘2rr 2 A ’ A ’ 

= { cos + ^*) cos M* + r ) + V’j]} 


fc = l /=! 

"CO 


rco 

a 2 / cos wrj(w) c/u> 

J -CO 

/ oo 

costur^o(^’) <Ej 

"00 


(4.5) 


It is easily seen that 4.5 is identical to the auto-correlation function Ro{t) of f 0 {t) 
because of the Wiener Khintchinc relationship: 

r<X> , . 

/? 0 (r) = / <f>o(uj)e MT du} (4.6) 

J -OQ 

Rq{t) — i cos u>Tcf>o(u)duj (4.7) 

J -CO 

where i is the imaginary unit. 
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1 his indicates that the whenever the ensemble average is considered, the simulated 
piocess f(t) possesses the auto-correlation, R(t) and the mean square spectral density 
4>(u) which are identical, respectively, with the target- auto-correlation Rq{t) and 
target spectral density <f> 0 (to). 

hor the present model the variables corresponding to time t, process /(£), and 
temporal frequency to, are the corresponding spatial variables namely, distance s, 
track height from mean h(t), and spatial frequency ft respectively. 

llie normalisation of the PSD (equation 4.3) gives y(ft) to be a gaussian distri- 
bution. 


/ 

tf(ft) = I 


n 2 



if |ft| < ft 0 


(4.8) 


0 otherwise 

The gaussian deviate can be generated by the transformation method [24] from a 
uniform deviate between 0 and 1. 


4.1.1 Generation of The Uniform Deviate 

The uniform deviate can be generated by many different methods. In this work two 
methods have been investigated, the subtractive method suggested by Donald E. 
Knuth [24], and the Data Encryption Standard (DES) based method [24]. The DES 
based method of generating pseudo-random numbers is preferred as the randomness 
of the result is much better than any other known algorithm, however this algorithm 
is much slower than the others. 

4.2 Prediction of The Track Profile 

The prediction of the track profile is done in the time-domain. The random process 
in the space-domain is sampled at discrete time steps. The distance covered by the 
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vehicle during a time interval is calculated using the space-time relationship 3.25. 
Ihe value of the track height at this distance is got by interpolation. 

1 he linear predictor algorithm [31] used is robust and can vary it’s characteristics 
with time to adapt to the changing input PSD. \ he predictor is a auto regressive 
filter defined as 

v 

X n = Yl Q k X n~ki (4.9) 

*=1 

and it’s prediction error is given as 

Cn = ?n ~ ^n* (4.10) 

Here a;* is the predicted value at instant n, x y| _* is the value at instant n — fc, and 
<*k are the weighting oooHicients which have' to be adjusted in order to give a good 
prediction. The predictive coefficients are chosen to minimise 

'X" 

E 'l 

n = 0 

the square error. If the source spectrum is not constant, it becomes necessary to 
re-optimise the coefficients continually. 


4.2.1 Optimisation of The Prediction Coefficients 

The total squared prediction error is 

" = £4 = £(*.-f>*.-0 J . 


(4.11) 




lb minimise M by choice of coefficients aj, equation Til has been differentiated with 
respect to the coefficients aj, 

p \ 


so, 


dM __ 

ddj 

n 

(■ 

p 



E «* E - r 

n ~jZ n ~k ^ 

£ 

k= 1 » 


n 
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(4.12) 


(4.13) 


\ 




It can be seen that fox' a doubly infinite sum, 

oo 

x n-jXn-k v Ar-f-1 * — ^ ^ X n Xn^j~^k 

n=— oo n— — oo 

so now equation 4.12 transforms to 



Rq&i + R\&2 “f* -^2^3 "f* * ' * ^ jRi 

i?iai + Ro a 2 + R1CI3 + • * * = i?2 

i?2 a l + Rl<^2 + Rq&3 + * * * = i?3 


(4.14) 


/here 

Rm /* ' %nZn+mi ni = 0 } 1, 2, • • • , p. 
n 

The above system of equation has to be solved for a k . In the present study it 
ias been solved using the Durbin and Levinson method [34] which requires much less 
;omputational effort than other methods. 

Since it’s impractical and undesirable to compute the infinite sums in the above 
equation a windowing procedure is used. The window is described as a set of weights 
spread over time as follows 

x' n = w n x n . 

Here w n is zero outside a finite range of interest. Now R m can be redefined as 

•' h+N-l 

R t= £ «+m’ m = M,2, ***,?. 

n=A 

The prediction co-efficients a k are recalculated when the mean square error in 
prediction becomes considerable large.. This results in fine tuning of the prediction 
coefficients with progress in time. 
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Chapter 5 


Presentation of Results 


The results of the numerical simulation of the model and the numerical values of the 
system parameters used for this simulation are presented in this chapter. The study 
was carried out for three different parameters namely landing sink velocity, landing 
glide velocity, and runway mean profile. 


5.1 The Model Parameters 

The values used are presented in table 5.1. These refer to a typical light fighter/trainer 

aircraft. 

The actual damping and spring relations of the aircraft were non-linear. They 
were approximated by two straight lines which formed a crude envelope of the original 

curve. 


5.2 Track Parameters 

The track is modelled as a sum of zero mean random process and a variable mean 
deterministic process. The parameters used for generating these different processes 


Symbol Description 

Value 

Units 

Ci 

Damping coefficient for shock absorber 

1.021x10 s 

Ns/m 

c 2 

Damping coefficient for tyre 

200.0 

Ns/m 

AT 

Spring coefficient for shock absorber 

1.0 xlO 6 

N/m 

AT 

Spring coefficient for tyre 

1.8xlO b 

N/m 

My 

Sprung Mass 

4133.33 

Kg 

EH 

Unsprung Mass 

77.0 

I<g 


Landing mass of aircraft 

8600.00 

Kg 

7 

Landing velocity 

75.56 - 77.0 

m/s 

Vs ink 

Sink velocity 

1.0 - 3.0 

m/s 

s 

Surface area of the wing 

38.4 

m 2 

D 

Diameter of the wheel 

0.66 

m 

I 

Moment of Inertia of tire wheel 

0.56 

I<g m 2 


Angle of attack at approach 

13.5 

deg 

|j 

Set angle of attack during ground roll 

6.0 

deg 

imi 

Stalling angle 

20.0 

deg 

tv* 

Angle of attack at touch-down 

16.85 

deg 

Eg 





Coefficient of friction drag at touch down 

0.0614 


Bfli 

Period of free roll 

0.038 

s 


Table 5.1: Model Parameters 
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arc as follows : 


e Zero mean random process : 

The zero mean random process was generated using the Shinozuka filter. The 
parameters used were 

1. Correlation constant for roughness, a = 0.005. 

2. Standard deviation of the track, cr = 0.01. 

3. Cutoff Nyquist frequency, 0; = 0.006. 

4. Number of simulation terms for the Shinozuka filter, N = 200. 

9 Variable mean deterministic process : 

The following four different mean profiles for the runway were generated : 

1. Flat runway : 

This is equivalent to a zero mean process alone. 

2. Inclined runway : 

This runway was assumed to have an inclination of about 1 in 1000, cor- 
responding to an angle of inclination of 0.06° I 

3. Stepped runway : I 

The step was assumed to be a AR-19 aluminium repair mat with a step! 
height of 0.038 m and a length of 16.45 m, placed at a distance of 100.0 m 
from the point of touchdown. 

4. Sinusoidal runway : ' i 

The sinusoidal runway was generated for a sine wave with a amplitude of 
0.05 m and a wavelength of 15.23 m. 

The study has been carried out for the cases listed in table 5.2 ! 
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Case 

Glide 

m/s 

Sink 

m/s 

Mean Profile 

Figure Numbers 

1 

75.56 

1.0 

Flat 

5.1, 5.9, 5.17, 5.25, 5.33, 5.41, 5.49, 5.57 

2 

75.56 

2.0 

Flat 

5.2, 5.10, 5.18, 5.26, 5.34, 5.42, 5.50, 5.58 

3 

75.56 

3.0 

Flat 

5.3, 5.11, 5.19, 5.27, 5.35, 5.43, 5.51, 5.59 

4 

76.32 

1.0 

Flat 

5.4, 5.12, 5.20, 5.28, 5.36, 5.44, 5.52, 5.60 

5 



Flat 

15.5, 5.13, 5.21, 5.29, 5.37, 5.45, 5.53, 5.61 

6 


1.0 

Inclined 

5.6, 5.14, 5.22, 5.30, 5.38, 5.46, 5.54, 5.62 



1.0 

Sinusoidal 

! 5.7, 5.15, 5.23, 5.31, 5.39, 5.47, 5.55, 5.63 

B 

75.56 

1.0 

Stepped 

5.8, 5.16, 5.24, 5.32, 5.40, 5.48, 5.56, 5.64 


Table 5.2: List of trial runs 


5.3 Presentation of Results and Discussion 

The results are presented in figures 5.1-5.64. Points are plotted 70 milli seconds 
apart, whereas the programme generates output after every 0.5 milli second. This 
results in considerable loss of information in the presentation, specially in the early 
stages where the effect of the landing impact is present. In latter stages where the 
system has stabilised the information loss is neglegible. 

The results for the following state variables are presented. 

1. Ground Reaction, 

2. Retarding Friction on the tyre, 

3. Sprung mass vertical displacement, 

4. Sprung mass vertical velocity, 

5. Sprung mass vertical acceleration, ' 

6. Unsprung mass vertical displacement, 

7. Unsprung mass vertical acceleration, 
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Zi 

^2 

Skid 
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cm 
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■ 


0.00 



0.00 

0.00 

-100.00 

0.00 

0.00 


0.70 


25.88 

51.18 

28.90 

70.65 

-66.34 

-12541.90 

-1.15 

■SjjSE 

-1.77 

0.197 

IgTtl 


7.68 

98.98 

43.55 

65.92 

5.91 

2599.37 



-1.81 

■mi 

yi flap 



143.50 

86.15 

61.29 

526.88 


3.29 

-10.43 


■mi 

■ElESl 



184.83 1 


56.77 


68498.30 


-10.58 

-1.80 

0.174 

preEBi 

32.45 

16.89 

223.02 



251.81 


-0.45 


-1.72 

0.201 

■fWl 

30.65 

0.00 

258.15 

52.34 

48.00 

giWEBl 

. 2635.34 

0.97 

-10.35 



mm 

■laatfl 






-10888.40 

0.53 


-1.71 

larazii 

■MMI 

mm 

6.32 

319.37 

47.50 

39.42 

165.71 

-846.36 


BEB3 

-2.35 


■asci 

29.97 




34.99 

-28.41 

8674.32 

mga 

roia 

-2.16 

0.173 

7.00 

8.12 

0.00 

368.36 

13.21 

30.64 




mm 



7.70 

16.99 

1.04 

388.29 

26.56 

26.34 

■EEI 


■RPI 



mil 

CO 

O 

28.41 

8.70 

405.24 

46.10 

22.07 

73.16 

2307.95 

0.26 

-10.93 

mam 

0.264 

9.10 

9.98 

0.49 

419.19 

17.90 

17.86 

58.91 

-33557.30 

■Em 




mail 

39.97 

19.17 

430.26 

66.52 

13.79 



■Em 




Bppini 

44.95 

13.86 

438.48 

68.01 

9.71 

16.48 

33807.50 

0.43 

-10.86 


0.239 

MWBftl 

26.08 

22.19 

443.95 

51.30 

6.22 

30.76 

11329.90 

-0.16 

-10.87 

-2.19 

■mi 

H 

0.00 

0.00 

447.28 

70.26 

3.13 

-175.86 


-0.21 


-2.19 

mm 

11.97 

26.70 

8.12 

447,49 

44.61 

2.78 

152.22 

-3869.63 

0.28 

-10.89 


0.922 


Table 5.3: Output for Landing Run on a Flat Runway, 
aircraft glide velocity = 75.56 ~ aircraft sink velocity = 1.0 y 


8. Braking Torque on the wheel, 

9. Skid of the wheel, and 

10. the forward ground -velocity of the aircraft. 

The graphs are plotted as functions of time. Distance travelled by the aircraft is 
also plotted as an alternative scale. However this is a non-linear scale as the forward 
velocity changes with time. 

A sample output is presented in table 5.3. This output is for the case of the 
aircraft landing on a flat runway with a glide velocity of 75.56 y, and sink velocity 
of 1.0 'f. 
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5.3.1 Ground Reaction 

Both the actual and the predicted ground reactions are plotted ( figures 5.1, 5.2, 5.3, 
5.4, 5.5, 5.6, 5.7, 5.8 ). It is observed that for all the cases the ground reaction is 
oscillatory in nature. A sharp peak in the initial response is present due to the landing 
impact. This peak is seen to increase in magnitude with increase in sink velocity, as 
is expected. The effect of the glide velocity ( figures 5.1, 5.4, 5.5 ) is neglegible on 
the initial impact response of the vehicle. However in the latter stages the reaction 
stabilises much faster with increase in glide velocity. The road holding capability of 
the aircraft is observed to increase with increase in both sink ( figures 5.1, 5.2, 5.3 ) 
and glide velocities in- the latter stages. 

The variation in the mean profile also affects the ground reaction. When compared 
to a flat track ( figure 5.1 ) the road holding is much poorer for a sinusoidal runway 
( figure 5.7 ), where the aircraft leaves the track after almost every peak. However 
this is much improved in the case of inclined ( figure 5.6 ) and the stepped runway ( 
figure 5.8 ). This could be due to the faster dissipation of energy in the earlier stages. 

The predicted reaction generated through an auto regressive filter, as discussed in • 
§4.2, is found to match the actual reaction quite closely. However sometimes at peaks 
it is unable to follow the actual reaction exactly. The prediction is seen to get better 
as time progresses due to finer tuning of the prediction coefficients and stabilisation 
of the response. 

5.3.2 Friction Force 

The friction force ( figures 5.9, 5.10, 5.11, 5.12, 5.13, 5.14, 5.15, 5.16 ) between the 
wheel and the ground is a function of fi the coefficient of friction and the ground 
reaction. The coefficient of friction itself is a function of the aircraft velocity and the 
skid of the wheel. 
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In the initial stages of the landing run the skid is maintained within a small 
range near the optimal skid value, corresponding to the maximum possible friction 
coefficient. The friction force therefore displays characteristics similar to the ground 
reaction ( figures 5.1, 5.2, 5.3, 5.4, 5.5, 5.6, 5.7, 5.8 ). In the latter stages, however 
the skid experiences stronger oscillatory changes and consequently the friction force 
is not able to stabilise as in the case of the ground reaction. This unsettled behaviour 
is prominent at the glide velocity 75.56 m/s in lower range of the sink velocity ( fig- 
ures 5.9, 5.10, 5.11 ). For the case of increased glide velocity ( figures 5.9, 5.12, 
5.13 ) and non-zero mean track profiles ( figures 5.14, 5.15, 5.16 ) the friction force 
characteristics are closer to the ground reaction characteristics. 

5.3.3 Sprung Mass Displacement 

The sprung mass displacement ( figures 5.17, 5.18, 5.19, 5.20, 5.21, 5.22, 5.23, 5.24 ) 
stabilises after an initial peak, due to the landing impact and goes to it’s static dis- 
placement value. The magnitude of the peak is observed to be dependent on the sink 
velocity as expected. Some kinks arc visible in the graph, which are manifestations • 
of the track profile undulations. 

For the case of sinusoidal track unevenness ( figure 5.23 ) the displacement ampli- 
tude is found to attain a maximum at 7.84 seconds after touch-down, corresponding 
to a velocity of 25.75 m/s, signifying resonance. The ground input frequency at 
this instant is found to be equal to the sprung mass frequency. After this point of 
resonance the displacement is seen to stabilise and go down to it’s static value. 

In the case of the stepped profile ( figure 5.24 ), the effect of the step is manifest 
in the form of a corresponding step in the sprung mass displacement. 
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5.3.4 Sprung Mass Velocity 

The sprung mass velocity ( figures 5.25, 5.26, 5.27, 5.28, 5.29, 5.30, 5.31, 5.32 ) shows 
exactly similar nature as the sprung mass displacement ( figures 5.17, 5.18, 5.19, 5.20, 
5.21, 5.22, 5.23, 5.24 ). 

5.3.5 Sprung Mass Acceleration 

The sprung mass acceleration ( figures 5.33, 5.34, 5.35, 5.36, 5.37, 5.38, 5.39, 5.40 ) 
shows clearly the effect of landing impact in the form of a peak with magnitude 
depending on the sink velocity. The variation in acceleration is seen to die down with 
the progress of the landing run. The effect of resonance in the case of the sinusoidal 
track profile ( figure 5.39 ) is not visible, however the presence of the step in the 
stepped runway ( figure 5.40 ) is clearly reflected by the acceleration graph. The 
effect of the inclined runway ( figure 5.38 ) is only to damp the acceleration sooner. 

5.3.6 Unsprung Mass Displacement 

The unsprung mass characteristics ( figures 5.41, 5.42, 5.43, 5.44, 5.45, 5.46, 5.47, 
5.48 ) are observed to be similar to that of the sprung mass. The sprung mass 
displacement is plotted along with the track height so that road holding is easily 
visible. Road holding is seen to be poor in the case of low sink velocities ( figures 5.41, 
5.42, 5.43 ) and is very bad for the sinusoidal runway ( figure 5.47 ). However for 
the sinusoidal runway the road holding improves after resonance has taken place 
in the sprung mass. The effect of resonance is also visible in the unsprung mass 
displacements as during this period the sprung mass exerts more effective weight on 
the unsprung mass than normal. 
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5.3.7 Unsprung Mass Acceleration 


The unsprung mass acceleration ( figures 5.49, 5.50, 5.5!, 5.52, 5.53, 5.54, 5.55, 5.56 ) 
does not show the effect of landing impact clearly. This could be due to the lower 
resolution of the graph. Variation in acceleration is seen to die down much faster with 
increase in glide velocity ( figures 5.49, 5.52, 5.53 ) and variable mean ( figures 5.54, 
5.55, 5.56 ) track profile. The effect of resonance for the sinusoidal track profile 
( figure 5.55 ) is also very small. The presence of the step in the profile ( figure 5.56 ) 
is found to be not easily discernible. 

5.3.8 Braking Torque and Skid 

The braking torque ( figures 5.57, 5.5S, 5.59, 5.60, 5.61, 5.62, 5.63, 5.64 ) is a function 
of the ground reaction and /z the coefficient of friction. This torque is the control 
variable that is being varied so as to maintain the skid in the optimal range of 0.16 
to 0.18. The brakes start acting only after the period of free roll. 

At higher values of the ground velocity a'small deviation in skid from the optimal 
range produces a large deviation in the coefficient of friction, whereas for lower values 
of ground velocity this deviation is much lower. Thus the range for optimal skid 
increases until at near zero velocity it extends up to full skid. It is seen therefore 
that initially the variation in skid ( figures 5.57, 5.58, 5.59, 5.60, 5.61, 5.62, 5.63, 
5.64 ) is small but as the ground velocity reduces the skid tends to go to higher values 
approaching full skid. 

Control on the brakes is applied keeping the above variation in the friction coef- 
ficient in mind. The braking torque is varied so as to keep the angular deceleration 
of the wheel a constant, if as a consequence of this the wheel goes out of the optimal 
skid range a correction is made to the value of the braking torque. This results in a 
wide variation of the braking torque with time. The minimum value of the braking 
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torque is obviously zero corresponding to the free rolling condition. The wheel is 
also allowed to roll freely when the aircraft leaves the track or balloons. 

Duration of continuous braking increases with increase in glide velocity ( fig- 
ures 5.57, 5.60, 5.61 ), specially in the latter stages of the landing run. Same is valid 
for the inclined track profile also ( figure 5.62 ). The stepped track profile ( fig- 
ure 5.64 ) has braking characteristics similar to the flat track ( figure 5.57 ). For the 
sinusoidal track ( figure 5.63 ), however, the braking torque is discontinuous and put 
to zero more often as the aircraft leaves the track more often. This results in poorer 
braking, even though the magnitude of the braking torque is much higher than in the 
other cases. 

5.3.9 Ground Velocity 

The ground velocity ( figures 5.1, 5.2, 5.3, 5.4, 5.5, 5.6, 5.7, 5.8 ) is seen to be almost 
linear except for the sinusoidal track ( figure 5.7 ) where it is rippled. Each ripple 
corresponding to the period when the aircraft leaves the runway. The deceleration in 
aircraft ground velocity at the start of the run is seen to be higher for higher values 
of sink velocity. This is due to the increase in retarding friction force which is a 
consequence of the increased reaction due to the landing impact. 

5.3.10 Braking Distance 

The expected ground roll for the aircraft used in this study is 740 m during landing. 
In the present study the ground roll ( figures 5.1 — 5.8 ) is found to vary from 416.6 
metres for the stepped track ( figure 5.8 ) to 493.0 metres for the inclined track 
( figure 5.6 ). This shows a marked reduction in the ground roll and the efficacy of 
the proposed braking system. The increased ground roll for the inclined runway is 
unexpected and is perhaps due to the weight component of the aircraft being neglected 
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in the present formulation. 

The ground roll is found to vary with the track mean profile ( figures 5.1, 5.6, 
5.7, 5.8 ) and is dependent very slightly on the sink ( figures 5.1, 5.2, 5.3 ) and glide 
( figures 5.1, 5.4, 5.5 ) velocities. It reduces marginally with increase in sink velocity' 
and shows a small increase with the increase in glide velocity. 
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5.1: Predicted and actual ground reaction for flat runway 
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Time ( seconds ) 

Figure 5.4: Predicted and actual ground reaction for flat runway , 
glide velocity = 76.32 f , sink velocity =1.0® 
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Figure 5.5: Predicted and actual ground reaction for flat runway 
glide velocity = 77.00 f , sink velocity = 1.0 — 
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Figure 5.6: Predicted and actual ground reaction for inclined runway 
glide velocity = 75.56 ~ , sink velocity = 1.0 — 
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Figure 5,7: Predicted and actual ground reaction for sinusoidal runway 
glide velocity = 75.56 j , sink velocity = 1.0 ~ 
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Figure 5.S: Predicted and actual ground reaction for stepped runway 
glide velocity = 75.56 ^ , sink velocity = 1.0 y 
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‘igure 5.9: Predicted and actual retarding friction for flat runway 
;lide velocity = 75.56 — . sink velocity = 10® 
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Figure 5.10: Predicted and actual retarding friction for flat runway 
glide velocity = 75.56 y , sink velocity = 2.0 y 
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Figure 5.21: Sprung mass displacement for flat runway 
glide velocity = 77.00 — , sink velocity = 1.0 — 
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'igure 5‘,26: Sprung mass velocity for flat 
,lide velocity = 75.56 —■ , sink velocity = 
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dgure 5.30, Sprung mass velocity for inclined r 
lide velocity = 75.56 ~ , sink velocity = 1.0 ~ 
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igure 5.43: Unsprung mass displacement for flat runway 
tide velocity = 75.56 f , sink velocity = 3.0 a 
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Figure 5.4S: Unsprung mass displacement for stepped runway 
glide velocity = 75.56 j- , sink velocity = 1.0 ~ 
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Figure 5.53: Unsprung mass acceleration for flat runway , 
glide velocity = 77.00 — , sink velocity — 1.0 — 
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Figure 5.53: Braking torque and skid for flat runway 
glide velocity = 75.56 > sink \elocit\ — *-»0 s 
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Figure 5.59: Braking torque and skid for flat runway 
glide velocity = 75.56 f , sink velocity = 3.0 f 





PIMS 


CM O CO <0 't CM O 

i4 T-t 



( ) anbJOi 6 uxm bj Q 

103 



PT»S 


04 

t-i 


o 

■H 


CO (0 c\j o 



( uiN>i ) 3 nt:) J01 SufM^JG 
104 


Time ( seconds ) 

Figure 5.62: Braking torque and skid for inclined runway 
glide velocity = 75.56 ~ 5 sink velocity = 1.0 — 
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Chapter 6 


Conclusions and Scope for Future Work 


Some salient conclusions gleaned from the present stud}' are presented in this chapter. 
The proposed scheme is very simple and a lot of work can be done by way of extension 
and otherwise to achieve better results. Some suggestions are made regarding the 
scope for such work. 

6.1 Conclusions 

The following conclusions are drawn from the results obtained from the various runs. 

• Marked reduction in the landing run can be obtained by using the present 
system of anti-skid braking. 

• The optimality of the proposed scheme is evident from the fact that the landing 
■ run is somewhat independent of both the sink and the glide velocities of the 

aircraft. 

• Any dominant frequencies in the runway profile adversely effect the landing 
run of the aircraft, therefore care should be taken to remove such dominant 
frequencies from prepared runway surfaces. 



• The introduction of a step early in the landing run, reduces the ground run. 
Concrete runways with differential settlement of slabs could therefore prove 
beneficial for the ground roll. 

• Higher landing impact gives higher initial deceleration, therefore the lift should 
be spoiled as early as possible in the ground run. 

• Inclined runways are not necessarily beneficial from point of view of reducing 
ground roll. 

6.2 Scope for Future Work 

The present work presents a very simple anti-skid braking system. Further work can 
be done by considering a more generalised system which accurately mimics the actual 
behaviour of the tyre and the landing gear. Systems with more degrees of freedom 
and with non-linear damping gears would be much nearer to the physical braking 
system. 

Active landing gears which optimise the suspension parameters for achieving 
higher and more steady ground reactions would definitely improve the optimality 
of this scheme by far. 

The track profile could be generated by using mathematically more advanced and 
complex methods e.g., chaos theory, for multivariate or univariate track processes . 
Further, simulation of the tyre’s rolling effect, which has been neglected in this study, 
can also be pursued. 

Other prediction methods which directly predict the state of the system, like the 
Kalman filter, instead of predicting the input as in the present study, would enhance 
the efficiency of the scheme. 


108 



Experimental verification of this scheme, which might produce some different re- 
sults would be desirable. 
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